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ABSTRACT 

We present a new approach for calculating the physical properties of highly ionized X-ray flows in active 
galactic nuclei (AGN). Our method relies on a detailed treatment of the structure, dynamics and spectrum of 
the gas. A quantitative comparison of our model predictions with the 900ksec Chandra/HETG X-ray spectrum 
of NGC 3783 shows that: (1) The highly ionized outflow is driven by thermal pressure gradients and radiation 
pressure force is less important. (2) A full featured dynamical model that provides a very good fit to the 
high resolution X-ray spectrum requires a multi-phased flow with a density power spectrum reminiscent of the 
interstellar medium. (3) Adiabatic cooling is an important factor and so is an additional heating source that 
may be related to the apparent multi-phase and turbulent nature of the flow. (4) The base of the flow is ~ 1 pc 
from the central object, in agreement with some, but not all previous estimates. (5) The mass loss rate is in the 
range 0.01 — 0.1 M yr _1 which is smaller than previous estimates and of the same order of the mass accretion 



rate in this object. 

Subject headings: acceleration of particles — ISM: jets and outflows - 
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1. INTRODUCTION 

Highly ionized gas (HIG) in active galactic nuclei (AGN), 
also known as "warm absorber", was first detected in the early 
1980s (Halpern 1984) and was observationally and theoreti- 
cally investigated in hundreds of papers (see e.g., Netzer et 
al. 2003, hereafter N03, and references therein). This com- 
ponent is believed to be photoionized by the central X-ray 
continuum source and is observed in some 50% of all type-I 
(broad line) AGN (e.g., George et al. 1998, 2000, Porquet, 
Reeves, O'Brien, & Brinkmann 2004). The poor resolution 
of the pie-Chandra & XMM instruments did not allow full 
investigation of the HIG properties and its possible relation 
to the UV outflowing gas (e.g., Crenshaw et al. 1999). The 
launch of Chandra & XMM allowed an in-depth analysis of 
the spectral properties of this gas component (e.g., Kaastra et 
al. 2000, Kaspi et al. 2000, Branduardi-Raymont et al. 2001, 
Lee et al. 2001, Steenbrugge et al. 2005). X-ray absorp- 
tion lines were detected in several sources and found to be 
blueshifted by a few x 100 kms~' relative to the systemic ve- 
locity (though much higher velocities have also been reported, 
e.g., Chartas et al. 2002). 

Detailed photoionization modeling of several AGN show 
outflows with a stratified density and temperature structure. 
Such models are partly phenomenological since they do not 
self-consistently account for the flow dynamics. Specifically, 
such models cannot be used to predict the velocity and deduce 
dynamically related quantities such as the mass loss rate. This 
information is crucial for understanding the effect of AGN 
flows on their environment (e.g., King 2003, Scannapieco 
& Oh 2004). A more complete approach requires full, self- 
consistent dynamical, photoionization, and spectral modeling 
of the flows. Today there are only a few qualitative works 
of this type (e.g., Arav, Li, & Begelman 1994, Murray et al. 
1995, Chelouche & Netzer 2003a,b). 
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This paper introduces a novel approach for investigating the 
properties of highly ionized flows. It utilizes the entire spec- 
tral information and employs state-of-the-art photoionization, 
dynamical, and spectral calculations. We apply this method to 
the X-ray flow in NGC 3783 and compare our model predic- 
tions to the 900 ksec Chandra/HETG spectrum of this source 
(e.g., Kaspi et al. 2002). The paper is organized as follows: In 
§2 we summarize the properties of NGC 3783 and its outflow. 
In §3 we use general arguments to distinguish between several 
plausible dynamical models. Section 4 outlines the new for- 
malism and §5 presents the self-consistent solution and shows 
a detailed comparison with the X-ray spectrum of NGC 3783. 
We elaborate on the implications of our results in §6 and sum- 
marize the work in §7. 

2. THE HIGHLY IONIZED GAS IN NGC 3783 

NGC 3783 is a Seyfert 1 galaxy that has been studied, ex- 
tensively, in almost all spectral bands. It has been the focus 
of a large UV and X-ray campaign with several already pub- 
lished papers (Kaspi et al. 2001, 2002, Behar et al. 2003, 
Gabel et al. 2003a,b, Krongold et al. 2003, N03). While 
certainly the best studied AGN in the X-ray band, NGC 3783 
is by no means unique and its low resolution ASCA spectrum 
resembles that of other AGN exhibiting warm absorption fea- 
tures (e.g., George et al. 1998). 

The bolometric luminosity of NGC 3783, L tot , can be esti- 
mated from its 2-10keV X-ray flux reported by Kaspi et al. 
(2002) and the spectral energy distribution (SED) presented 
in N03. This results in L tot ~ 3 x 10 44 ergs -1 which we as- 
sume to be the long-term average luminosity of the source. 
The mass of the black hole, Mbh, has been estimated from re- 
verberation mapping to be ~ (3 ± 1) x 10 7 M and the size of 
the broad line region (BLR) ~ 10 16 cm (Peterson et al. 2004; 
see however Kaspi et al. 2000 for lower values). Assuming 
solar composition for the gas we define the Eddington ratio, 
r = 1.18L to t/iEdd. (e.g., Chelouche & Netzer 2001) where 
^Edd. is me Eddington luminosity. For NGC 3783 T ~ 0. 1 . 

The superb 900 ksec Chandra/HETG X-ray spectrum of 
NGC 3783 allowed the first in-depth study of its HIG outflow. 
The results were thoroughly discussed in Kaspi et al. (2001, 
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2002), Krongold et al. (2003) and N03. Here we consider the 
Kaspi et al. (2002) and N03 results as a basis for our model. 
The main results of these papers can be summarized as fol- 
lows: The HIG is outflowing with v ~ 1000 km s and the 
absorption line profiles exhibit extended blue wings. The ve- 
locity (either centroid or dispersion) does not show significant 
correlation with ionization level. The HIG consists of at least 
three distinct ionization components ([low, intermediate, and 
high]) whose ionization parameters, U ox (defined as the ratio 
of the photon density in the energy range 0.54-10keV to the 
hydrogen number density) are log(f/ OT ) = [—2.4, — 1.2, — 
0.6], respectively. Each component has a different column 
density, N H , given by, log(N H ) = [21.9, 22.0, 22.3] for the 
above three values of U ox . The N03 analysis shows that none 
of the HIG components responded to the large, observed con- 
tinuum variations. Recombination time arguments suggest, 
therefore, a lower limit on their distance of [3.2,0.63,0.18] pc, 
respectively. The line of sight covering factor of the flow, Ci os , 
is between 0.7 and 1 .0. The global (Ait) covering factor of the 
flow, C, obtained from the emission lines is ~ 20%. The mass 
loss rate, M, was estimated to be ~ 75 eC M yr _1 , where e 
is the volume filling factor. For e = 1, it exceeds the expected 
mass accretion rate by roughly two orders of magnitude. 

3. GENERAL MODEL CONSIDERATIONS 

Several models have been proposed to qualitatively explain 
the acceleration of highly ionized gas in AGN (e.g., Krolik 
& Begelman 1986, Balsara & Krolik 1993, Konigl & Kartje 
1995, Bottorff, Korista, & Shlosman 2001, Chelouche & Net- 
zer 2001). The models we consider here assume that the flow 
is driven either by radiation pressure force, thermal pressure 
gradients, or a combination of the two. The equation of mo- 
tion in this case is 



dv _ 1 
dr p 



n e a T Uoi ( M 1 \ dP 
Anr 2 c V T J dr 



(1) 



where r is the distance from the ionizing source, v the veloc- 
ity, p (n e ) the gas density (electron number density), Oj the 
Thomson cross-section and P the thermal pressure (P = pv 2 
where v s = ^2k B T/fim H ~ 166(T /IQ 6 K) 1 / 2 kms" 1 and 
fXtriH is the average mass per particle; pi ~ 0.6 for fully ion- 
ized solar metallicity gas). M is the force multiplier defined 
as the ratio of the total radiation pressure force to that due 
to Compton scattering (see Chelouche & Netzer 2003a and 
references therein). For the densities associated with the HIG 
(<C 10 12 cm -3 ), the ionization structure of optically thin flows 
depends, to a good approximation, only on U ox . In this limit, 
M is a function of U ox and the optical depth in the line per unit 
thermal velocity width (e.g., Chelouche & Netzer 2003a). 

The modeling scheme considered here combines detailed 
photoionization and spectral calculations (some aspects of 
which were discussed by N03) with a simplified treatment of 
the gas dynamics. As shown below, the unique set of obser- 
vational constraints allows us to quantitatively investigate, for 
the first time, which of the aforementioned dynamical models 
is most relevant to the HIG outflow in NGC 3783. 

3.1. Radiation pressure driven flows 

Non-rotating radiation pressure dominated outflows require 
M > F-\ Fi gure^ suggests that this condition is difficult to 
meet for solar metallicity, optically thin HIG, assuming the 
known F and SED. Specifically, this condition is not met for 
the intermediate and high ionization components presented in 
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FIG. 1 . — The force multiplier, M, as a function of U ox for the SED used 
by N03. Optically thin and thick (where M is an average over the cloud) 
solar composition gas is shown in solid lines. The gravitational term for 
NGC 3783 is shown in horizontal dash-dotted line. The three vertical lines 
mark the HIG components of N03. Note that unless the flow is highly metal 
rich (top dashed curve), radiation pressure force is smaller than gravity for 
two of the three N03 solutions. 



N03, where we find that M < T -1 for U ox > 5 x 10~ 3 . There 
is evidence that BLR metallicities can exceed solar (e.g., 
Hamann & Ferland 1993, see however Shemmer et al. 2004). 
In our case, the metallicity should exceed about 30 x solar for 
radiation pressure force to overcome gravity for U ox < 0.3. 
Such metalicities have never been observed in AGN. 

Radiation pressure force can exceed gravity close to the 
foot-point of rotating flows. As shown by Chelouche & Net- 
zer (2001), HIG flows must be launched within or just outside 
the broad line region (~ 10 16 cm) in order to reach the ob- 
served velocities. Unless the line broadening in NGC 3783 is 
unrelated to the gas dynamics, such models cannot account 
for the width of the line profiles given the much larger dis- 
tances implied by the Behar et al. (2003) and N03 analyses. 
The above considerations do not take into account the effect 
of dust on the flow dynamics (e.g., Scoville & Norman 1995, 
Everett 2002) since none has been detected in NGC 3783. We 
conclude that radiation pressure force is unlikely to drive the 
flow in NGC 3783 to the observed velocities. 

We note that HIG flows may not be steady-state phenom- 
ena. In this case, the HIG may have been accelerated in the 
past (when M was larger or the object brighter) and we ob- 
serve it now at its coasting or even decelerating phase. Al- 
though plausible (e.g., Chelouche & Netzer 2001), such a 
model adds little to our understanding of the general phe- 
nomenon. 

3.2. Thermal pressure driven flows 

Acceleration by thermal pressure gradients (see Parker 
1958) has been suggested to drive extremely ionized high 
temperature gas in AGN (e.g., Begelman, McKee, & Shields 
1983, Krolik & Begelman 1986, Balsara & Krolik 1993 
and Woods et al. 1996). In particular, Balsara & Krolik 
(1993) and Woods et al. (1996) presented the results of de- 
tailed, time-dependent hydrodynamic calculations of thermal 
pressure driven flows for a wide range of initial conditions. 
The conclusion is that highly ionized gas near its Compton 
temperature may be driven to terminal velocities of up to 
~ 2000 km s . Here we investigate the possibility that such 
a mechanism drives the ionized gas in NGC 3783. 
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FIG. 2. — Normalized velocity profiles of thermal pressure driven poly- 
tropic flows. r c is the location of the sonic point (see text). Larger y (i.e., 
cases where adiabatic cooling is more important) result in lower velocity for 
the same sonic point. Also shown is the density profile for a y = 1 isothermal 
flow (dashed line) which resembles the hydrostatic density profile for r < r c . 
The column density per unit velocity for an isothermal flow at a fixed ion- 
ization parameter (°= r(v)~ 2 dr/dv; dash-dotted curve) illustrates the way the 
opacity decreases with velocity (see §4.3) 



For a spherically expanding flow, 

poerV 1 . (2) 

In the limit of M <C V (see figure [T), and assuming a poly- 
tropic flow with index y (y = 1 for isothermal gas) we can 
rewrite equation^as 



1 dv 
v dr 



1 



yvj 



2yv 2 



n e a T L tot 
Anpr 2 cT 



(3) 



Super-sonic steady-state winds that are launched sub- 
sonically should cross the sonic point at r c = GMbh/2/v 2 , 
where v(r c ) = \/yv s {r). At this point both the numerator and 
the denominator vanish thus the equation of motion is defined 
at all radii. Noting that Mbh x ^Edd. x ^tot/r we find 
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10 6 K 
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(4) 



For a given y, wind solutions are self-similar and scale with 
v s and r c . Several such solutions for different values of y are 
shown in figure [2] The terminal velocity of such flows is of 
the order of (or a few times larger than) the sound-speed at the 
critical point (i.e., the escape velocity). For isothermal flows 
(7=1) the terminal velocity is a few times the sound-speed 
(~ 4v s at r ~ 10 2 r c ) while for polytropic flows with 1 < y < 
5 /3 the terminal velocity is lower due to adiabatic losses (e.g., 
Lamers & Cassinelli 1999). Polytropic flows where y < 1 
imply excessive heating and result in higher flow velocities. 

Also shown in figure[2]is the density profile for an isother- 
mal flow. For r < r c the density profile is similar to the hydro- 
static case. The density beyond the sonic point declines more 
rapidly compared to the hydrostatic case. The column density 
of the flow, from some r to infinity, is 



N H (r) = / dr'n H {r'). 



(5) 



For super-sonic isothermal flows, v <* ^/ln(r) and Nh ~ 
rn H (r). 

Returning to NGC 3783, we notice that the temperature of 
the high ionization component found by N03, combined with 



the mass of the central object in this source and assuming 
7=1 gives (equation |4} r c ~ 3pc. This distance is consis- 
tent with the three lower limits on the distance found by N03. 
Using equation[5]and the density of the hot component found 
by N03, we get Nh ~ 10 22 cm~ 2 which is similar to the col- 
umn density of this component obtained by N03. Thus, the 
conditions in this source are consistent with a picture where 
the highest ionization gas is in a state of thermal pressure 
driven flow. As shown below, there are other observational 
constraints that are satisfied by such a model. 

3.3. The density structure of the flow 

The continuity condition (equation[2j implies U ox v. This 
has been shown to result in high ionization absorption line 
profiles that are more blueshifted and broader compared to 
low ionization lines (Chelouche & Netzer 2002). Such trends 
are not observed in NGC 3783 (Kaspi et al. 2002). In 
fact, lines of different ionization states are remarkably similar 
which indicates similar dynamics. It is unlikely that different 
ionization components are launched with different initial con- 
ditions and end up showing the same line profiles. A more 
likely explanation is that all flow components are co-spatial 
and are accelerated together to their observed velocities (a 
similar conclusion has been reached by N03 based on the ap- 
parent pressure equilibrium between the different phases). 

The co-existence of several gas phases is a well known phe- 
nomenon in the interstellar medium (ISM) and in molecular 
clouds. A multi-phase gas naturally occurs for a turbulent 
medium where the relation between the density p and the 
physical scale B, follows a powerlaw distribution, p=<^ (see 
§4 for the full definition). Such dependences have been ob- 
served on large (1 — 10 3 pc; e.g., Elmegreen & Scalo 2004) as 
well as small (10 15 — 10 18 cm; Deshpande 2000) scales. 

Consider a small section of the flow over which both r and 
v are constant. Since Nh(^) « p^ x and U ox « p _1 , we 
get N H « t/<£ where a = -Q3 + l)/]8 (cf. Krolik & Kriss 
2001). Thus, there is a simple local relationship between 
the column density and the ionization parameter involving the 
size distribution parameter of the medium. This dependence 
holds also for geometrically thick flows and the integrated col- 
umn density provided the self-similar scaling holds for all r 
(see §4.3). 

Returning to NGC 3783, we find that the N03 three com- 
ponent solution implies j3 ~ —0.8. Such a value of j3 usu- 
ally characterizes turbulent media like the ISM and molec- 
ular clouds and is consistent with recent simulations (e.g., 
Boldyrev, Nordlund, & Padoan 2002). 

Our model assumes the following range of densities at ev- 
ery location r 



P=Pn 



(6) 



where p m ; n and ^ max are the (location dependent) minimum 
density and maximum length, respectively, and j3 is a constant 
to be determined later. We require that the surface of a sphere 
with radius r is fully covered by all phases in the density range 

[PminjPmax]? 

Jdt(p) 



x2 = Anr 2 . 



(7) 



This results in 
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FIG. 3. — The density profile of a multi-phase flow with a fixed density 
contrast. Also shown are lines of constant ionization parameter. Clearly, the 
expansion of an optically thin (to continuum absorption) flow results in an 
increase of its mean level of ionization. 
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FIG. 4. — The mass loss rate, M, as a function of the densi ty powerlaw index 
j3 for different values of the density contrast (eauation ll II . M increases with 
decreasing ft and increasing density contrast since the fraction of the denser 
regions of the flow is larger in those cases (see text). 



where ^ max is defined over the surface of the sphere. Assum- 
ing that the density contrast, p ma x/Pmin, is a constant of mo- 
tion then <^ max oc r . 

The mass loss rate due to all flow components can be cal- 
culated by defining the mean density over a (thin) spherical 
shell, 



(p(r)> = 4^r^ ;r 



(9) 



which may be expressed, using the above relations, as 



(P(r)) = 



27T(j3+2) V r 



6. 
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n \ GB+2)/j3' 
Pmax \ 



The mass loss rate is then given by 



(10) 



M = A-nr 2 (p(r))v = 



8ftr 2 p min V 1 - (Pmax/Pmin) 



(j3+2)/j3 



/3+2 1 - (Pmax/Pmin) 2 ^ 



(ID 

This conforms with the mass loss rate of a single phase spher- 
ical flow in the limits J3 — * 0, p max = p m ; n . Assuming a diver- 
gence free flow, mass conservation implies 



Pmin'" v = Const. 



(12) 



This continuity condition is identical in form to that of a 
single-phase medium. It means that the highest ionization pa- 
rameter, f/™ ax , corresponding to the lowest density p m i n , is 
proportional to the velocity. By construction there are fixed 
density and ionization parameter contrasts at every location 
(Pmin/Pmax = U™ n /U^ ax ). An example of a density profile 
of a multi-phase flow is shown in figure [5] together with sev- 
eral iso-ionization parameter curves. Clearly, the mean level 
of ionization increases as the gas expands. 

For a given velocity profile, the mass loss rate depends on 
the density contrast and p\ For a fixed p m i n , smaller j3 and/or 
larger density contrasts result in denser regions having larger 
scales and a larger contribution to M. The differences are 
small for /3 > — 1 since the mass fraction of the denser regions 
is negligible (figure^. 

The effect of inhomogeneity on the flow dynamics can be 
approximated by noting that all components in the N03 solu- 
tion are in a rough pressure equilibrium. Under such condi- 
tions one must solve for the combined dynamics of all phases 



since their coupling is strong and the momentum is quickly 
(relative to the dynamical timescale) distributed throughout 
the flow. The dynamical effect can therefore be included by 
considering their contribution to the inertial mass. The mean 
(over volume) density is ~ 3/Q3 + 3 ) p m i n - In this approxima- 
tion, the critical point is given by 



r e = [(l+j8/3)y]- 



GM BH 

2v? : 



v(r c ) = v,V(l+/V3)y. (13) 



In the above analysis we have neglected the effect of drag 
forces between the different phases. This is justified if 
such condensations are constantly forming and evaporating 
on short timescales compared to the dynamical time and the 
mean net effect is averaged-out. Thus, the flow considered 
here is very different from other works in which discrete, 
dynamically independent clouds, rather than time-dependent 
condensations, are assumed (e.g., Chelouche & Netzer 2001, 
Everett 2002). 

3.4. Thermal balance in a multi-phase expanding medium 

As demonstrated in several works (e.g., Begelman, McKee, 
& Shields 1983, Balsara & Krolik 1993, Woods et al. 1997), 
adiabatic cooling can have an important effect on the thermal 
equilibrium and, hence, on the dynamics of continuously ex- 
panding Compton-heated winds in low luminosity AGN. Here 
we investigate the effect of such cooling on the thermal struc- 
ture of the multi-phased HIG in NGC 3783 where the temper- 
ature is below the Compton temperature and the medium may 
be turbulent. 

Neglecting the small contribution due to velocity gradients, 
the adiabatic cooling term for a spherically expanding contin- 
uous flow is 

2 v 

Aadiabatic - 2pv f -. (14) 

r 

The N03 model of NGC 3783 did not include adiabatic cool- 
ing. Considering the observed properties of NGC 3783 (r = 
10 19 cm and v ~ 1000 km s _1 ) and including the process in 
the ION thermal calculations, we find that adiabatic cooling 
lowers the maximum temperature by more than an order of 
magnitude (see figure [5}- Since we do not have a numeri- 
cal hydrodynamic code, we have tested this result against the 
numerical calculations of Woods et al. (1996) under similar 
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FIG. 5. — The effect of adiabatic cooling on the stability curve of pho- 
toionized gas in NGC 3783. In all models the gas is assumed to flow with a 
velocity of 1000 km s , the distance is 10 19 cm, and the SED and luminos- 
ity (3 x 10 44 erg s _1 ) are those defined in N03 (see §2). The thick solid line 
is the N03 stability curve with their assumed, r = 10" cm and no adiabatic 
cooling. The thin solid line has the same parameters but includes adiabatic 
cooling. Additional stability curves that include adiabatic cooling were cal- 
culated by varying one parameter at a time: r, L. The parameter that was 
changed is denoted next to the curves. 



conditions (such as the Eddington ratio, SED, and gas veloc- 
ity; see their figure 14) and find them to be in good agreement. 
Thus, the numerical approach adopted here gives the correct 
temperature and expansion velocity which is required by our 
calculations. The effect is smaller when the ionizing flux is 
higher since heating is proportional to L tot /r 2 (see the indi- 
vidual curves in figure[5j. The ionization structure of the flow 
and, hence M, are less affected by adiabatic cooling. 

Adiabatic cooling also alters the stability curve for the gas. 
The three N03 solutions still lie on the stable branches of the 
stability curve yet the pressure of the high ionization compo- 
nent is smaller by roughly a factor 3 so that strict pressure 
equilibrium between all components is not maintained. 

The flow considered here is multi-phased and turbulent (see 
§6.3). In such a case there is an additional internal energy 
source entrained in the large eddies. This energy is trans- 
ferred to smaller and smaller scales until it dissipates (e.g., via 
Coulomb collisions) and serves as an additional heat source 
for the gas. The corresponding heating rate is 

,3 



'turbulence : 



v tmb 

r 



(15) 



where v t m-b is the turbulence velocity (see Bottorff & Ferland 
2002). Typically, v tur b is of the order of v s and the flow veloc- 
ity, v (e.g., Shu 1992). The size of the largest eddy is com- 
parable to the flow extent and, hence, to r. Hence, heating by 
turbulent dissipation can offset the adiabatic cooling. Below 
we consider several models where the two terms play different 
roles. 

4. A THERMAL PRESSURE DRIVEN FLOW MODEL FOR NGC 3783 

In §3.2 we have shown a qualitative agreement between 
the predictions of a multi-phase, thermal pressure driven flow 
model and the observations of NGC 3783. Section 3.3 out- 
lined the physical model in detail. Here we carry out a more 
quantitative investigation of the model and apply it to the out- 
flowing gas in NGC 3783. The main aim is to give detailed 
spectral predictions and compare them to the high resolution 



X-ray observations. We use the simplified, steady-state treat- 
ment of the flow dynamics presented in §3.2 and a sophisti- 
cated scheme for the calculation of the transmitted spectrum 
through the flow. As explained, the changes in the gas tem- 
perature and the attained velocities are in general agreement 
with more sophisticated hydro-calculations. 

Here the flow is assumed to evaporate from some mass 
reservoir (e.g., the accretion disk or the putative torus), ex- 
pands as it is accelerated to cover a fair fraction of the AGN 
sky, and crosses our line-of-sight at some finite radius (to be 
determined later). Such flow configurations are often assumed 
in ID modeling of AGN flows (e.g., Murray et al. 1995) and 
have been qualitatively confirmed by detailed hydrodynami- 
cal simulations of energy driven flows (e.g., Balsara & Rrolik 
1993; see also Woods et al. 1996). 

4.1. Calculation Scheme 

We begin the calculation by specifying the initial condi- 
tions. These include the location, ro, where the flow enters 
our line-of-sight (and need not be equal to r c ), t/™ ax (ro), and 
t/™"(ro) which also determine the density contrast through- 
out the flow. An additional parameter is the density powerlaw 
index, /3 . The other parameters, L t ot and T, are fixed at their 
observed values. 

The next step involves photoionization and thermal calcula- 
tions of non-LTE gas in statistical equilibrium (this is justified 
since the dynamical timescales are larger than the recombi- 
nation timescales by several orders of magnitude; Rrolik & 
Kriss 2001, N03). These calculations are carried out using 
ION 2004, the 2004 version of the ION photoionization code 
(e.g., N03, with the addition of including new di-electronic re- 
combination rates for iron; see Netzer 2004) which includes 
all important heating and cooling mechanisms and is suitable 
for exploring the physics of HIG in photoionization equilib- 
rium. The code includes also the effect of adiabatic cooling 
on the thermal structure of the highest ionization phase (see 
equation Il4li . Once the temperature profile of the high ion- 
ization phase is calculated, we approximate it by a polytropic 
relation (T °< ) which provides a good approximation for 
the cases considered here, and solve the equation of motion. 
(We note that, for the relevant parameter space, radiation pres- 
sure force acting on the low ionization phases of the gas has 
a negligible effect on the global flow dynamics). The process 
iterates several times until v(r) and / converge. 

The next step is to calculate the ionization and thermal 
structure of all flow phases at location r. This is done for a 
flow section of length 8r in the radial direction. 8r is taken as 
the scale over which the flow velocity changes by 10km s . 
This ensures high resolution spectral predictions for the ab- 
sorption lines. For this purpose we use a discrete version of 
the continuous density distribution (equation[6} that includes 
50 phases corresponding to roughly 0.1 dex intervals in U ox . 
When calculating the transmitted spectrum through sections 
of the flow we assume that equation[8]holds also for the line- 
of-sight direction albeit with a different normalization (spher- 
ical condensations or "clouds" are not the general case under 
such conditions; e.g., Blumenthal & Mathews 1979). This 
means that 

p = p min (l70^ (16) 
where is the length scale in the radial direction and 

Pmax 



C/ 
s max 



l - 



(17) 
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FIG. 6. — The effect of j3 on the total column density distribution as a 
function of ionization level (i.e., U 0K ). The linear part of the curve is sensitive 
to p (see text). Steeper density profiles (smaller /3) result in the flow being 
dominated by low ionization gas. 
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FIG. 7. — The effect of j8 on the transmitted spectrum. Steeper density 
profiles (smaller j8) result in low ionization gas contributing more to the gas 
opacity with a stronger suppression of the soft X-ray flux. 



B, 1 is determined by the somewhat arbitrarily defined Sr and 
is therefore independent of | . The ratio | /S, 1 is constant, by 
construction, for all phases, and is their geometrical aspect ra- 
tio (e.g., Blumenthal & Mathews 1979). This arbitrariness in 
the shape and the number of condensations (see the appendix) 
has no bearing on the results presented in this work which rely 
on observationally deduced quantities, and is related to our 
ignorance concerning the hydrodynamical properties of such 
flows and the exact shape (e.g., the number of absorption sys- 
tems) of the X-ray line profiles. 

The contribution of all phases to the line and continuum 
opacity at each section of the flow is calculated assuming full 
line-of-sight coverage of the continuum source (see below). 
The transmitted continuum then serves as the ionizing spec- 
trum for the photoionization calculation of the successive sec- 
tion at r + 8r. The velocity shifts at each point are taken into 
account when calculating the transmitted spectrum. The cal- 
culation continues to large r and produces a high resolution, 
theoretical X-ray spectrum. 

The global covering fraction of a phase with density p in 
the range [p m m(r),p ma x(r)} at a distance r is given by 

C(p, 



4nr 2 t, l {p) r VPmin(r) 



(18) 



Clearly, the covering factor, expressed in this way, can be 
much smaller than unity for the denser phases and can sig- 
nificantly effect the transmitted continuum through the flow. 
However, this is only a geometrical definition and the more 
relevant quantity is the velocity dependent covering factor, 
C(p , v), which takes into account the gas motion and thus en- 
ables the proper treatment of the opacity (e.g., Gabel et al. 
2003). We note that, in a turbulent flow, there is no one-to- 
one correspondence between v and r. Turbulence broadens 
the lines and increases C(p,v) (see also §5). 

4.2. Dependence on initial conditions 

The solution of the self-consistent model depends on sev- 
eral initial and boundary conditions. The dependence on the 
initial conditions of the problem, j3 , t/™ ax , t/™ m , and ro is dis- 
cussed here for the case of isothermal flows. Qualitatively 
similar results are obtained also for polytropic flows. 

We first explore the dependence of the column density dis- 
tribution on U 0K and j3. As discussed in §3.3 and shown in 
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FIG. 8. — Absorption line profiles as a function of t(vq) (upper-left), 
I'turb (upper-right), vq (lower left), and v, (lower-right; here, i»o = v., is as- 
sumed). The model marked by a blue line in all panels corresponds to 
t(v ) = 10 2 , Vturb = 1: vo = 200km s _1 , and v., = 200km s _1 . All line pro- 
files exhibit an extended blue wing. 



figure |6j affects the column density distribution such that 
smaller j3 results in larger columns of low ionization gas. The 
expansion of the flow and the imposed boundary conditions 
result in rapid decline of the column density for extreme val- 
ues of U ox . Different column densities result in different trans- 
mitted spectra. This is shown in figure where we find that 
smaller j3 result in more opaque flows. The effect can be dra- 
matic even for 0.2 change in j3 in cases of large integrated 
columns. 

Smaller density contrasts, i.e., larger t/™ n at a fixed U™ x , 
result in smaller columns of less ionized gas (figure [6} and 
therefore lower opacity and higher transmitted soft X-ray flux 
(i.e., softer spectra, figure^. 

The distance, ro, where the flow crosses our line-of-sight 
with velocity vo, affects the total column density of the flow. 
The effect is manifested as a change in the X-ray "color". For 
example, smaller ro imply larger densities, more opaque gas, 
and harder X-ray spectrum. This is especially important if our 
line-of-sight crosses the subsonic section of the flow (ro < r c ) 
given the nearly hydrostatic structure of the gas. 

We next consider the effect of initial conditions on the ab- 
sorption line profiles and assume, for demonstrative purposes, 
that lines are produced throughout the flow by a single ioniza- 
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tion phase (i.e., over a narrow range in U 0K ). In such a case, 
the optical depth as a function of velocity (for a phase with a 
given t/ ox ) is 



t(v) oc / dv'r(v') 



- 2 {dr(v') 



V dv' 



exp 



(19) 



where we have assumed a Gaussian kernel for line broaden- 
ing due to e.g., a turbulent velocity field. The dependence of 
the line profile on t(vo) (the optical depth at vo), v tur b, and v s 
is shown in figure [8] The calculations show that all line pro- 
files exhibit extended blue wings due to the flow expansion 
with the most extreme line asymmetries occurring for opti- 
cally thin lines (e.g., f(vo) ~ 1). Larger v tur b broadens the 
lines and makes them more symmetrical with respect to vo. 
vo naturally causes a line shift such that larger vo will result 
in more blueshifted, narrower lines for the same v(r). When 
calculating the effect of the sound-speed (of the highest ion- 
ization component), v s , on the line profile we have assumed 
vo = v y which is relevant to super-sonic flows that are opti- 
cally thin in the continuum. As shown, larger v v imply larger 
terminal velocities and considerably broader lines. Obviously, 
the observed line profiles may have a more symmetrical shape 
depending on the spectral resolution of the instrument. 

5. FITTING THE 900 KSEC CHANDRA SPECTRUM OF NGC 3783 

We resort to an iterative scheme to fit the observed high 
resolution X-ray spectrum of NGC 3783. We first try to con- 
strain the flow dynamics. This is accomplished by taking the 
theoretical line profiles, normalizing the column densities to 
agree with those obtained by N03, and convolving the calcu- 
lated profile with the known instrumental profile (Kaspi et al. 
2002). A comparison between the predicted and the observed 
line profiles defines a narrow range of r c (hence, v s and con- 
versely C/™ ax X v o and Vt ur b. An initial estimate for /3 comes 
from the scaling of the column density with ionization pa- 
rameter (see §3.3). This allows us to employ the next step in 
our fitting scheme which is a more systematic and time con- 
suming search in a localized region of the parameter space, 
involving complex spectral calculations (section 4.2). 

The calculated line and continuum spectra take into ac- 
count the probability for our line-of-sight to intersect any 
given phase at every location in the flow (see §4.2). We as- 
sume that the finite thermal and turbulent broadening results 
in C(p,v) = 1 (the opacity at every velocity bin is properly 
scaled so that the integrated column density is conserved). 
The best-fit model is obtained by comparing the detailed theo- 
retical spectrum to the observed spectral shape and to individ- 
ual absorption line profiles. All calculations assume a photon 
continuum slope, r = 1.6, which best characterizes the hard 
X-ray spectrum of this source. 

Considering only radiative heating and cooling, and adia- 
batic cooling, we find a set of models (hereafter adiabatic 
models) which are consistent with the spectral shape and 
whose critical point is about 0.5 pc from the center and 
ro = r c . The density contrast of the flows spans more than 
four orders of magnitude at every location (e.g., t/™ m ~ 
10~ 4 , t/™ ax ~ 0.4 at r ~ 1 pc). Such flows reach a veloc- 
ity of about 600 km s and the absorption line profiles, con- 
volved with the instrumental resolution, are shown in figure[9] 
Clearly, the predicted lines are narrower and less blueshifted 
than observed. A better fit is obtained by requiring that the 
flow crosses our line of sight at large distance (~ 1 .5 pc) and 
the line profiles are broadened by a turbulent velocity field. 




FIG. 9. — Comparison of the calculated Si XIV A6.182 line profile with 
the third order MEG data of NGC 3783. The solid line shows an adiabatic 
model which crosses our line-of-sight at its sonic point at 0.5 pc. Such a 
model cannot account for the observed blueshift or width of the line. A non- 
turbulent flow that crosses our line-of-sight at 1.5 pc has a more blueshifted 
centroid but a small equivalent width (dashed line). A similar but turbulent 
( v turb=l50 tmr 1 ) A ow P rov ides a better fit to the data (dash-dotted line). 



For Si XIV A6.182, a v tul -b ~ 150 km s _1 provides a better fit 
to the data (figure[9). 

The observational indication for the existence of turbulence 
suggests that the thermal equilibrium of the gas may be af- 
fected by it. Due to the uncertainties associated with the 
combined effects of adiabatic cooling and turbulent dissipa- 
tion heating we study a case where A a( n a batic = ^turbulence and 
radiative processes are the only important heating and cooling 
mechanisms. In this case, the high ionization, volume filling 
phase of the flow may be considered isothermal (e.g., Begel- 
man, McKee, & Shields 1983). We refer to it as the isother- 
mal model but note that all components which are not fully 
ionized have radial temperature dependence. 

The best-fit isothermal model is shown in figure ma (a 
similar fit is obtained for models using either of the full 
line-of-sight coverage methods described in §4.2). Clearly, 
the model is acceptable and is qualitatively similar to the 
best-fit model of N03. Inspection of individual absorption 
line profiles (figure [TT1 shows qualitative agreement with re- 
spect to the line centroids and widths. The derived mass 
loss rate is M ~ 10- (I8±0 - 6 C M yr" 1 [for C ~ 0.2 (N03), 
M r-j 0.01 — 0.1 M Q yr _1 ], i.e., considerably smaller than the 
N03 and the Behar et al. (2003) estimates. The model also 
gives detailed predictions for the column densities of various 
ions. These are given in table 1 for the more abundant ele- 
ments. Clearly, the stratified structure of the flow results in 
a wide range of ionization levels having comparable column 
densities. We discuss the implications of these results in §6.2. 

The range of parameters characterizing acceptable models 
are given in table 2. The quoted uncertainties refer to differ- 
ences between the adiabatic and isothermal models, the mass 
of the central object, and the degeneracy in the spectral fit. 
Specifically, the error in the value of each parameter repre- 
sents the maximum deviation from the best-fit value for which 
a good fit can be obtained by varying all the other parameters 
in the table. Thus, uncertainties in different parameters are 
not independent. 

6. DISCUSSION 



8 



Chelouche D. & Netzer H. 




5 5.5 6 6.5 7 7.5 8 8.5 




c 9 9.5 10 10.5 11 11.5 12 12.5 13 

o 

o 




13 13.5 14 14.5 15 15.5 16 16.5 17 




17 17.5 18 18.5 19 19.5 20 20.5 

Rest wavelength (A) 



FIG. 10. — Full fit to the 9()0ksec Chandra/HETG X-ray spectrum of NGC3783. This is a pure absorption model and line and continuum emission is not 
included. The model provides an acceptable fit over most bands but does not reproduce well the shape of the UTA near 15-16A. The solid line assumes full 
line-of-sight coverage of the source and the dashed line in the bottom panel an 80% coverage. The mismatch near 5.5A is due to a well known calibration effect 
(N03). 



TABLE 1 
Ion Column Densities (log(cm~ 2 )) 



Element 


I 


II 


III 


IV 


V 


VI 


VII 


VIII 


IX 


X 


XI 


XII 


XIII 


XIV 


H 


16.5 




























He 


13.7 


19.3 


























C 


7.4 


12.2 


15.3 


16.6 


17.8 


18.0 


















N 


6.5 


11.5 


14.9 


16.5 


16.8 


17.4 


17.5 
















O 


6.8 


11.9 


15.5 


17.0 


17.4 


17.5 


18.1 


18.2 














Ne 


6.0 


11.3 


14.8 


16.2 


17.0 


17.2 


16.8 


16.8 


17.4 


17.6 










Mg 


7.8 


10.4 


13.6 


14.7 


15.3 


16.1 


16.4 


16.6 


16.7 


16.7 


17.1 


17.2 






Si 


6.5 


10.2 


12.2 


13.4 


14.9 


15.8 


16.2 


16.5 


16.6 


16.7 


16.7 


16.6 


17.1 


17.2 



Column densities of various ions for our best-fit model. Boldface values mark deviations by 
more than a factor two from those deduced from observations and reported by N03. 



We have introduced a new approach for modeling highly X-ray flow in NGC 3783 is that of a multi-phased, geomet- 
ionized flows in AGN. The emerging physical picture for the rically thick flow that is thermally evaporated from a parsec 
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FIG. 11.— Best-fit model to the MEG and HEG 1st and the MEG 3rd order 
absorption line profiles (model parameters are given in table 2). The spectral 
model (solid line) was convolved with the instrumental resolution relevant to 
each case. Note the clear blue asymmetry of all predicted, non-convolved 
line profiles (dashed red-lines). 

TABLE 2 
Best-fit model parameters 



Parameter 


Value 


Critical distance r c [10 18 cm] 


2.0±1.5 


Line-of-sight crossing distance ro [10 I8 cm] 


5±2 


Maximum ionization parameter at ro [log^™ 3 ")] 


0.5 ±0.3 


Minimum ionization parameter at ro [log(J/™ ln )] 


-3.7± 1.5 


Density spectrum index (j3) 


-1.2±0.5 


Mass loss rate [Mq yr _1 ], 


2q-U.8±0.6£* 



For C ~ 0.2 (N03), M is in the range 0.01 - 0. 1 M a yr _ 



scale region. The density range (~ 10 2 - 10 6 cm- 3 at r ) oc- 
cupied by the flow is much larger than obtained by Rrongold 
et al. 2003 and N03 and implies a density of the warm ab- 
sorbing phase which is, on average, lower than that suggested 
by Reynolds & Fabian (1995) and Bottorff & Ferland (2001). 
Our model has extremely ionized dynamical components that 
are not detected spectroscopically. The least ionized compo- 
nents are also absent in previous works since it is difficult to 
disentangle the ionization-level-column-density degeneracy. 

6.1. Is the flow turbulent? 

The above analysis suggests that the outflowing HIG may 
be turbulent for several reasons: 1) The flow is multi-phased 
with different phases at a rough pressure equilibrium. 2) 
The density spectrum (i.e., /3) is similar to that found in the 
ISM and in molecular clouds which are known to be turbu- 
lent media. 3) The Reynolds number of the flow is Re ~ 10 7 . 
This is obtained by assuming non-magnetic plasma and tak- 
ing the scale of the largest eddies to be ro, the velocity v c , 
and the density of the highest ionization phase (see e.g., equa- 
tions 3.25 & 3.30 in Lang 1999). Such a medium is capable 
of sustaining turbulent motion. The density or range in our 
model spans roughly five orders of magnitude (correspond- 
ing to U™ x /U™ n ). This agrees with the naive theoretical ex- 
pectation that the maximum density is that at which energy 
cascade is terminated and the energy is dissipated. This oc- 
curs when the density range spans roughly ~ 31og(7?e)/4 ~ 5 
decades (e.g., Shu 1992); i.e., consistent with the results ob- 
tained from the spectral fit. 

Our model does not explain the onset of turbulence in the 
flow. One possibility is that it is triggered by the putative 



central star cluster within the central parsec of NGC 3783. 

6.2. The UV-X-ray connection 

Table 2 suggests that the X-ray flow in NGC 3783 may 
also produce strong low ionization UV absorption lines of 
H I X 1216, C IV A 1548, and O VI A 1035. Our predictions are 
consistent with the lower limit on the column density of C IV 
deduced by Crenshaw et al. (1999). The model suggests that 
the global covering factor is smaller for lower ionization gas 
phases (equation I18>. This, and the fact that the UV con- 
tinuum emitting source is likely to be larger than the X-ray 
source, imply that the UV line profiles may be prone to par- 
tial line-of-sight coverage effects. This may be manifested 
either as multi-component absorption (i.e., not all flow phases 
are represented at every velocity bin; see §4.2) or as non-black 
saturation (e.g., Gabel et al. 2002). The current resolution of 
the X-ray gratings and the signal-to-noise level of the obser- 
vation are not high enough to tell which of these scenarios is 
relevant to NGC 3783. 

The dynamical timescale of the flow, ro/v s , is about 
10 4 years. Gabel et al. (2003) reported a considerable veloc- 
ity shift (~ 50 km s~ 1 ) of the highest velocity component over 
less than a year. Clearly, the two timescales do not agree. It 
is however possible that a filament moving across our line-of- 
sight rather than a real deceleration of the flow is responsible 
for the observed effect (see also Gabel et al. 2003). The rel- 
evant timescale in this case is the time it takes a filament to 
cross the surface of the UV emitting region. Assuming this 
region to be of order 1QR S , where R s is the Schwartzschild 
radius (~ 10 13 cm in the case of NGC 3783) and the velocity 
across the line-of-sight to be of the order of the thermal speed, 
we get a crossing timescale of ~ 0.3 year, i.e., consistent with 
observations. We conclude that either the high velocity decel- 
erating UV component is not related to the X-ray flow or else 
the apparent deceleration is caused by line-of-sight effects. 

6.3. Extended emission 

Some Seyfert 1 galaxies show evidence of extended, diffuse 
X-ray emission (e.g., NGC 3516; George et al. 2002). This 
is also common in Seyfert 2 galaxies. A natural question is 
whether the extended HIG flow is partly responsible for such 
emission through the scattering of continuum photons into our 
line-of-sight. The expected scattered X-ray luminosity assum- 
ing pure electron scattering is 



T x ■ 

-^scat - 



z e CL t( 



(20) 



where % e is the Compton optical depth and L[ ot the total X-ray 
luminosity. It is possible to relate z e C to the total mass loss 
rate of the flow, M, via 



GjM 
AntriH Jr, 



r~ 2 v~ l dr, 



(21) 



where r e is the distance beyond which L* cat is calculated. This 
gives 



^~6x 10~ 6 



M 



M yr- 



kpc\ / 10 kms~ 



(22) 



The extended (r e ~ 0.5 kpc) X-ray flux from NGC 3783 is 
< 10 36 ergs -1 (George I., private communication), i.e., con- 
sistent with model predictions (see table 1). 
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6.4. Relation to other AGN components 

So far, we have not considered the source of the outflow- 
ing gas. Our model suggests that, for NGC 3783, it is located 
within the central ~ 1 pc. A possible mass reservoir on such 
scales is the putative molecular torus (e.g., Krolik & Begel- 
man 1988; Krolik & Kriss 1995; and Netzer et al. 2002). Tak- 
ing the torus parameters from Netzer et al. (2002), we obtain 
that the torus mass is ~ 10 5 M Q and that HIG flows may be 
sustained at their current mass loss rates for ~ 10 6 — 10 8 years. 
If, as suggested by Krolik & Begelman (1988), the torus is 
also responsible for feeding the AGN, than this timescale also 
corresponds to the present activity life-time. 

The narrow line region (NLR) lies at the outskirts of the 
HIG flow and hence is exposed to a modified type-I AGN 
continuum. Similar conclusions (albeit with a somewhat more 
diminished UV fluxes compared to the model presented here) 
have been reached by Alexander et al. (1999) based on the 
analysis of narrow emission lines in NGC 4151. A physical 
interaction between the HIG outflow and the NLR is unlikely 
to be energetically important given the small mass loss rate 
derived here (cf. Laor 1998). 

6.5. Model limitations and extensions 

Despite the success of our model in explaining the observed 
X-ray features there are several notable problems. 

6.5.1. flow dynamics 

In its simplest form (without turbulent line-broadening or 
dissipative heating) our model under-predicted the line widths 
and blueshifts. One possibility, discussed in §5 and §6.1, is 
that the flow is turbulent. Alternatively, higher flow velocities 
and better agreement with the data may be obtained if the flow 
is launched closer in to the central source. This requires, how- 
ever, that the mass of the central object be smaller by a factor 
of a few compared to that reported in Peterson et al. (2004; 
see also §6.5.2). Another possibility is that the flow dynamics 
reflects a higher mean past luminosity compared with the one 
deduced from the present observations. In both cases, adia- 
batic cooling is less important (see §3.4) and radiative accel- 
eration is more effective (since Y is larger). Our calculations 
show that in such cases the terminal velocity is slightly larger 
(by ~ 300 km s ) than that of the isothermal model (see §5) 
and the fit is qualitatively similar. 

6.5.2. flow location 

The value of ro derived here is smaller than the lower limit 
derived by N03 on the distance of the multi-phased flow but 
is larger than the upper limit for the highest ionization com- 
ponent by itself. Nevertheless, the model is consistent with 
the lack of appreciable equivalent width (EW) variations of 
the silicon and sulfur lines since the multi-phase flow is ge- 
ometrically thick. Note also that r c and ro depend linearly 
on the mass of the central object which is somewhat uncer- 
tain. For example, Kaspi et al. (2000) report M BH ~ 10 7 M , 
i.e., a factor three lower than the new Peterson et al. (2004) 
estimate. A model having similar spectral properties would 
have r c ~ 0.13pc and ro ~ 0.5 pc (compared with the present 
r c = 2pc and ro = 5pc) for such low Mbh (equations- 
Reeves et al. (2004) reported XMM-Newton Epic-pn ob- 
servations showing significant EW variations of Fe XXIII and 
Fe XXIV absorption lines near 2A. The lines are consis- 
tent with the highest ionization N03 model yet, according 
to Reeves et al, they suggest the gas lies within the central 



0.02 pc. We have checked this suggestion within the general 
N03 model. Assuming the reported line variability is real and 
using the mean N03 continuum luminosity, we find that the 
photoionization timescale of Fe XXIII gives an upper limit of 
~ 0.2 pc, basically identical to the value denoted in N03. This 
ro is, however, inconsistent with the lower-limits on the lo- 
cation of the two more neutral, higher density components. 
Thus, the Reeves et al. (2004) observation presents a real 
difficulty to our multiphase model where all components are 
present at all locations. It remains to be seen whether the high 
ionization gas component reported by Reeves et al. is part of 
the outflow in NGC 3783. 

6.5.3. thermal instability 

Our model assumes a continuous range of densities and 
therefore ionization levels and temperatures. However, some 
phases of our model lie on thermally unstable branches of the 
stability curve (see figure |5}. In this case the gas may un- 
dergo non-linear compression or expansion following some 
initial perturbation and end up in a thermally stable equilib- 
rium. This will effectively change the assumed p — | distribu- 
tion. The appropriate modeling of such conditions is beyond 
the scope of this paper. 

6.5.4. spectral deviations 

Our model predicts lower flux levels then observed at soft 
X-ray energies. This discrepancy is eliminated if one takes 
into account partial covering (N03). Such an effect is, in fact, 
expected since the low ionization gas is likely to have lower 
global coverage. A continuum leakage of 20% (80% flow 
coverage) improves our fit and correctly predicts the depth of 
the (saturated but non-black) O VIII lines. Another discrep- 
ancy is evident near 15 — 16A where our model predicts too 
much absorption by highly ionized iron lines. This may be re- 
lated to the di-electronic recombination rates of some of these 
ions (cf. Netzer 2004). Also, figure^|does not include emis- 
sion features from the flow which results in under prediction 
of the flux in several bands (e.g., near the La line of Ne IX at 
13.6A; see Kaspi et al. 2002). 

Kaspi et al. (2002) showed that some X-ray line pro- 
files exhibit two absorption components. Our model assumes 
a smooth flow and, therefore, smooth line profiles. Multi- 
component absorption may be accounted for in two ways: 1) 
It is possible that the high velocity component seen in Kaspi 
et al. (2002) is due to a flow which is launched closer in to 
the central source where the temperature of the highest ion- 
ization component is higher and the escape velocity larger. 2) 
Absorption deficit at some intermediate velocities (e.g., due 
to a "missing" phase at some location; see §6.2.2). Modeling 
such multi-velocity flows is beyond the scope of this paper. 

7. SUMMARY 

We presented a novel method for modeling the properties 
of highly ionized flows in AGN. We performed state-of-the- 
art photoionization, dynamical, and spectral calculations, and 
compared our results to the 900ksec Chandra/HETG spec- 
trum of NGC 3783. The model is consistent with several in- 
dependent observational constraints including the continuum 
and line spectrum, spectral variations (or the lack of), and the 
upper-limit on extended emission from this source. Our ma- 
jor conclusions can be summarized as follows: 1) The flow in 
this source is geometrically thick and is driven primarily by 
thermal pressure gradients. 2) The flow is multi-phased with 
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a density-size spectrum which resembles that of a turbulent 
medium. Adiabatic cooling and turbulent heating are likely to 
be important in the energy balance. 3) The gas is located be- 
yond the broad line region and within the narrow line region. 
Its origin may be related to the putative torus. 4) The mass 
loss rate is of order of the mass accretion rate in this object 
and the kinetic luminosity is very small (< 0.1%) compared 
to the bolometric luminosity of the source. 5) The same flow 
may also account for some of the UV lines observed in this 
source. 
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APPENDIX 

ALTERNATIVE DERIVATION OF THE FLOW EQUATIONS 

It is possible to derive the flow equations using the volume filling factor formalism adopted by Arav, Li, & Begelman (1994) in 
their modeling of outflowing BAL regions. The Arav et al. model assumes confined clouds and may not be suitable for describing 
the highly ionized gas condensations considered here. Nevertheless, a comparison between the two approaches is useful. 

Define a volume filling factor by 

E(P) - £max ( — ) (Al) 

We require that all flow phases fill the entire available volume, i.e., / de = 1, to obtain 

/ Mmax 1 

e, 



1 



(A2) 



For a divergence free spherically expanding flow, 



a 2 f Pmm , / \ 4nr 2 p min V 1 - (p m ax/Pmin) T,+1 , 

M = A%rvi dpe(p) = ; ; — —, —. — — — =const., (A3) 

Jp mm + l 1 - (Pmax/Pmm)" 
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which, for a flow with a fixed density contrast, results in 

Pmin^V = Const. (A4) 

This is identical to equation[21and suggests that the continuity condition is not due to our choice of formalism. 

Consider now a volume V filled with a number of different density phases where a phase of density p occupies e(p) of 
the volume. If that phase is comprised of n identical spherical clouds then the volume of each cloud is n~ 1 e(p)V. The total 

surface covered by such clouds is ~ n x (n- l eV)) 2/3 = n 1 / 3 (eV) 2/3 . Thus, the global (An) covering factor due to this phase is 

C,(p)oc„V3 £ (p)2/3. 



